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Abstract 

Perturbative corrections to b — > c transitions at zero recoil are considered in the 
two-step matching scheme. The matching coefficient for the b — > c currents from 
the intermediate effective theory (between the scales nib and m c ) to the low-energy 
effective theory (scales below m c ) has been found with two-loop accuracy. The 
next-to-next-to- leading logarithms has been summed in the leading order in m c /mb. 
Higher-order corrections are estimated in the large-/?o limit. 
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1 Introduction 



Precise measurement of V c b is one of important tasks in the investigation of 
the CKM mixing matrix [1]. Two methods are currently used: inclusite and 
exclusive. Here we shall discuss the exclusive method. 



Weak b — > c transitions are described by QCD currents cTb with T = 7 M and 
757 M . If we are only interested in their matrix elements for pb = rribV + k, 
p c = m c v' + k', where the residual momenta k, k' are small as compared to 
the heavy-quark masses m Cj b, then we may expand these QCD currents in 
HQET operators of the form <v • ■ • b v with the appropriate quantum numbers 
divided by powers of m Cj b- The coefficients are calculated by matching. This 
is descussed, e.g., in the review [2] or in the textbook [3]. In this paper we 
consider the zero-recoil case v' = v. Then, in the v rest frame, 
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cy b = r] V c v b v + 0{l/m 2 c b ) , 
c 7 57& = rj A c v ^b v + 0{l/m 2 c b ) . 



(1) 



There are no l/m c ^ corrections [4]. We use MS renormalization with anticom- 
muting 75; then the vector and axial currents in QCD do not depend on the 
normalization scale \x. The HQET currents with v' — v are also /^-independent, 
and hence the matching coefficients r)v,A do not depend on any normalization 
scale. 

In the exclusive method, the differential rate of B — > D*W is extrapolated 
to the zero-recoil point. Thus, F(l)|Vrf,| is measured, where F(l) = i]a[1 + 
0{l/m 2 cb )}. Matrix elements of higher- dimensional HQET operators in the 
power correction are estimated in several ways, giving about 4% error in 
F(l) [1]. This error can be reduced to 1% level by unquenched lattice sim- 
ulations currently in progress. Therefore, it is desirable to know r/ r with the 
accuracy better than 1%. 

These coefficients can be obtained by matching QCD matrix elements with 
those in HQET: 



Here mj )C are the on-shell heavy-quark masses. In this way, we know the exact 
dependence of r]^ on r = m c /mb; however, we cannot sum terms with leading 
powers of L = \og(mb/m c ) in all orders of perturbation theory. The two-loop 
correction r/^ has been calculated in [5] and confirmed by an independent 
re-calculation [6]. The three- loop correction has been recently calculated [7] 
in the equal-mass limit m b = m c . 

Terms with the highest powers of (3 in all orders of the perturbative series 
were summed in [8]. The leading infrared renormalon in zero-recoil matching 
coefficients is at u — 1. This means that high-order behaviour of perturbation 
series for r]r is much better than in, e.g., the pole mass (where it is at u — 1/2). 
We don't have to worry about it yet. 

In this single-step matching approach, the evolution of a s and the relevant 
operators between the scales mi, and m c is not taken into account properly. In 
this interval, there are 4 active flavours; b is heavy while c is still light, and the 
evolution of operators is governed by HQET (in particular, a^(fx ) is not a 
natural quantity). And with (m c ) / (m b ) 1.56, these evolution effects 
are of order 1, not small corrections. Therefore, it is natural to take them into 
account from the outset. 

This is achieved by performing the matching in two steps. First, we match 




(2) 



2 



QCD with HQET-1, where b is static and c is dynamic, at /x ~ mj. Second, 
we match HQET-1 with HQET-2, where both b and c are static, at /x ~ m c . 
In this way, we obtain i] Vt A as series in r: 



r^r = C r (m b )G(m b ,m c )E(m c ) + 0(r) . 



(3) 



In practice, few terms of this expansion can be found, because the number of 
relevant operators in HQET-1 grows fast with the dimensionality. However, 
in each term of the expansion in r, we can sum leading, next-to- leading, etc., 
powers of L. The leading term in the r-expansion is known at the next-to- 
leading logarithmic order [9]. The 0(r) term is known at the leading logarith- 
mic order [10]. The 0(r 2 ) term is discussed in [11]. The two approaches can 
be combined: we subtract several first terms of expansion of (2) in r and add 
the RG-improved results instead. 

The QCD -> HQET-1 matching coefficients 



are known at the next-to-next-to-leading (NNL) order [12,13]. The running 
factor G(mb, m c ) is determined by the anomalous dimension 7 of the heavy- 
light current in HQET-1. Until recently, it was only known at two loops [14,15]. 
Now it is has been calculated at three loops [16], i.e., also at the NNL order, 
using the methods of [17]. Finally, the HQET-1 — > HQET-2 matching coeffi- 
cient 



is currently known at one loop [9] only. In the present paper, we shall find the 
two- loop coefficient e 2 , thus completing the NNL order calculation of t] V ,a at 



2 Summing NNL logarithms 

By requiring that the two-step matching result (3) (with the three-loop 7 [16]) 
reproduces the O(r ) term in r\y\ [5], we have two independent derivations of 




(4) 




(5) 



O(r ). 



the matching coefficient E(m c ): 



E(m c ) = 1 - 4CV 



47T 
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+ c, 



C F (-20vr 2 log 2 + 14C 3 + 73C 2 + ^) 

+ C A (l07r 2 log2 - 41C 3 - 6C 2 + ^) 

/ 233\ / 31 

+ Tp [-17(2 + ^ J + T F m (-4^ - ^ 



' ot s {m c ) 
An 



4 a s (m c ) , ( a s (m c ) 
1 - - V ' + 4.00 v ; 



3 vr 



7T 



(6) 



Here n\ = 3 is the number of flavours lighter than c; a^ 4 ) (/z) is written simply 
as a s (n). This result can also be obtained by a direct calculation; this will be 
discussed in a separate publication. 



Now we can sum the logarithms (up to NNL ones) to all orders: 
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i] r = x 



-7o/(2/3 ) 



An \ An 



(7) 



Here 



a s (m b ) ' 



(8) 



a s (/u) and (3 are for 4 flavours. The NL correction is 



n V,A 
Cp 



4 \ 1 21 



+ 



23 

8C2- T )CV-2(C 2 + 4)C A 



x — 1 



+ 



3C A (llC F + 7C A )(x- 1) 



(9) 



where the upper number in braces is for V and the lower one for A. The NNL 
correction is 



n V,A 

c F 



29\ 



3C2 + 7T7 ) x ~ 3 C 2 - 7T7 



1 53 



24/ 



+ 



24 I 205 
247. n 1.1 



/ 247 
^-207r 2 log2 + 3C 3 + — C2 + 



96 



x 



143 I _z 
231 I 16 



+ {^}^ a -{i}^-{r}^s{s 
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+ 



(l0vr 2 log2-|c3-fc 2 - 



41\ 2 



+ 



+ 



( 483 \ 
20( 4 + 9< 3 - 50(2 + -jg-J x 



^70 16 713 



+ 



„> 57 > 77 > 629s 



x- 1 



+ (6C 4 - 33C 3 + y C2 + ?f) (x + l)C\\ 

+ |^(32C 2 - 23) 2 (x - l)C% - (7( 2 + ^) (x + l)Ci 

17, 479\ Jrt> 193, 1 f 1281 
40C 4 + y C 2 + -), + 40 C4 + -C 2 -3 1545 



clc 
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5C 4 + 33C 2 - — J x - 5( 4 + C2 - { 



CfC\ 



x-l 



4 3C A (11C F + 7C A ) { (4(2 - y) (x - 1)C 2 F + \{x + l)C\ 



+ 



5\ , 27 



C f Ca 



x — 1 



+ 



9C F C1(11CV + 7C A ) 2 (x - if 



(10) 



Numerically, 



W,A = X 



6/25 



1 - 1.10588x4 



f 0.89412 ) \ a s {m b ) 
I 1.56079 J J vr 
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+ I 3.94156* 2 + I °- 98879 > , 



2.44780 ) \ f a s (m b ) \ 2 
1.72604 J X \ 8.22155 J ) \ n ) 



'ID 



3 Higher orders in the large-/3 limit 



It would be extremely difficult to sum the N 3 L logarithms. However, it is 
easy to find the (3$ term in ■ The highest powers of (3q in all orders of 
perturbation theory for 7 were found in [12]: 



7 = 



C F (3 



3 + 2/? 



3 fa B{2 + 0, 2 + - /3)r(3 + /3) 



= -3C, 



A) 



83 



(12) 



where /3 = /5 o; s /(47r). The matching coefficients Cr(m b ) were also considered 
in [12]: 



C VjA (m b ) = \-C F ^- 



+ 3C 2 



I 47/16 
I 445/48 

-1751 
7993 



P 2 + 



(13) 



Here we use the same method for E. The solution of the RG equation can be 
written as 



70/(2 A)) 



\o; s (m c ) J 

7(«s) 7o \ 



K(a s ) = exp y 



^2/3(a s ) 2(3 J a s 



(14) 



The RG invariant E is, at the first order in l//?o, 

1 \ 



E = 1 + 



2 / ' 



(15) 
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where (3 = /9 o; s (m c )/(47r). We obtain 

(5/3)u r»r(2 - 2u) i-u-u* 



S(u) = 3C F 



T(3 - u) 1 + 2u 



2u 



(16) 



Note that the leading infrared renormalon is at u = 1. Collecting all this 
together, we obtain 



E{m c ) = 1 + C>|- 
Po 



-4 + (3 C2 + H), 



+ 7C3-11C2 



7655 
432 



/3 2 + 



(17) 



This reproduces the (3q term at two loops (6), and gives us the (3$ term at 
three loops. 



Using this result, we obtain 



1,(3) 
a V,A 

Cp 



(6C3-HC2 

+ O((3 ) ■ 



3703 \ 3 

X ~ 6Cs 



13 U + J_ 751 

21 I ^ 216 1-4121 



/?o 2 
(18) 



Numerically, this means adding the term 

36.3 121 I \ ( a s (m b ) \' 



- 40.5461a; J + 



88.0131 



7T J 



inside the bracket in (11). Of course, this is an estimate only, because the 
terms O((3 ) and (9(1) are not known. 



4 What is known about the three-loop coefficients? 



Re-expressing a s (mb), a s (m c ) via af'(/io), we can write i]r in the form (2). 
From our NNL result, we obtain the L 3 , L 2 , and L terms in the three-loop 
coefficients r/^ . The (3q part of the non-logarithmic term can also be found. 
We can even make one more improvement, and obtain the (3% part of rj^ 
exactly as a function of r (including power corrections). The highest powers 
of (3q in all orders of perturbation theory for t]y were found in [8]. They have 
the form similar to (15), with (3 = /3oQ; s (/xo)/(47r) and 
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Sv( M )=6C^/')- r(u r ) . r ' 1 -. 2tt ) 1 -»- U ' [(1 - 2u)R v - flj , 

1 (3 — u) 1 + 2u 

r(3 — u) 1 + 2u 



x 



3(1 — 2u)(l — u — u )Ro - (3 - u + w 2 )i?! 



(19) 



where 



• i (l-2u)L 
Slllll 

i?o = cosh(L-u) , i?i = —j — 

sinh ^ 



(20) 



This gives us the (3% parts of as exact functions of r. 
Collecting all this together, we obtain 
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Vv'a 
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(21) 



Non-logarithmic terms of orders 0((3q) and 0(1) are only known [7] in the 
equal-mass limit r = 1: rffi' = and 



r,f /256 32, 4o 64 2l 2o 160 2 , 

440 184 448 2293 \ 
yC4 + ~~yC3 + ~~yC2 y- )^ F ^° 
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1696 9152 39296 16958\ 
+ — C 4 - -^Cs - ~^-(2 + -^-)CfT f 

( 9728 1216, , 1024 2l „ 17728 2l n i/ , ft/ . 
+ — a 4 - — log 4 2 + — 7r 2 log 2 2 + -^n 2 log 2 - 160Cs 
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+ — C2C3 - — C 4 - -^Cs ~ -gj— C2 + ~^-)C A T F 

+ (^Cs-^C 2 -^)Tl, (22) 

where a 4 = Li4(l/2). Of course, for the real value of r this is an estimate only. 
Assigning a ±100% theoretical error to it, we obtain a numerical estimate of 
the three-loop contribution to (2): 



1.81 ± 2.89^-W - 815 U 2 + M 27 61 ' 1 76 - 9 



73.0 
124.6 



l-rJ I 9.76 I \ I 45.6 | | 113.5 | 1 - r 

±3.9 



'«£ 6) (*>)' 



7T 



(23) 



(we used coth(L/2) = (1 + r)/(l - r)). 



5 Numerical results 

The program RunDec [18] was used for numerical estimates: 

m fe ^4.7GeV, m c «1.6GeV, r^0.34, L«1.08, 

a s (m b ) « 0.216, a; « 1.56, af 5 \fx ) « 0.258 . (24) 

With these numbers, the two-loop results [5] give 



rj v « 1 + 0.1899^- + 0.5416 — « 1-019, 

7r y 7r y 

rj A « 1 - 0.4768^ 1.3813 — ^ 0.952. (25) 

7T \ 7T / 



To combine the results of the two-loop calculation [5] with our resummation, 
we subtract from [5] its O(r ) term and add the resummed result (11) instead. 
This results in the corrections 

Ar] V « 0.021 , Ar] A ~ 0.029 , (26) 
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to be added to the results of [5]. If we also add the (3$ term (18) of the N 3 L 
logarithmic correction, these numbers change significantly: 

Ar] V « -0.049 , Ar] A « -0.059 . (27) 
It is more instructive to discuss the a 3 terms (23) in (2). They give 



A Vv n -(9.7 ±3.9) ( a ° 5) ^o) \ ~ _ .005± 0.002, 



«i 5) (/io) 



3 



At]a ~ —(19.6 ± 3.9) ~ -0.011 ± 0.002 . (28) 



7T 



The (3q terms in (21), i.e., the naive nonabelianization prescription [12], give 
the coefficients of (a s /7r) 3 in (28) equal to 0.9 and —14.0 for rjv and t]a] in the 
first case, even the sign is wrong. The authors of [7] estimated these coefficients 
as ± 3.9 and —11.1 ± 3.9, on the basis of their equal-mass calculation (our 
error estimate is the same as their one). Our results differ from their ones by 
about two estimated errors. The reason for this difference is the terms with 
powers of L which have been investigated in the present paper. 

Our final results 

rj V = 1-014 ± 0.002 , rj A = 0.941 ± 0.002 (29) 



have theoretical errors of about 0.2%; the central values differ from the previ- 
ous results [7] by about two estimated errors. 
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